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ABSTRACT 


During  the  1955  and  1957  Test  Operations  at  the  Nevada  Test  Site  (NTS),  masses  and 
velocities  were  determined  for  more  than  20,000  objects,  such  as  glass  fragments  from 
windows,  stones,  steel  fragments,  and  spheres,  which  were  energized  by  blast  winds  resulting 
from  nuclear  explosions.  Following  the  field  tests,  a  mathematical  model  was  devised  to  help 
explain  quantitatively  the  experimental  results.  This  model  required  certain  aerodynamic -drag 
information  in  regard  to  the  displaced  objects.  It  was  the  purpose  of  the  study  outlined  in  this 
report  to  determine  the  necessary  drag  properties  for  the  objects  by  means  of  drop  tests. 

In  addition  to  the  objects  mentioned  above,  small  laboratory  animals,  mice,  rats,  guinea 
pigs,  and  rabbits,  were  used  in  the  drop  tests.  The  data  obtained  from  these  te^ts  were  ex¬ 
trapolated  to  estimate  the  drag  properties  for  man,  and  the  results  compared  favorably  with 
data  from  other  sources.  Also  a  method  was  developed  to  estimate  the  average  drag  proper¬ 
ties  of  man  from  his  total  surface  area,  assuming  that  every  possible  orientation  of  a  straight, 
rigid  man  with  respect  to  the  wind  was  equally  likely. 
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NOMENCLATURE 
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a  =  acceleration 
Cd  =  drag  coefficient 

D  =  sphere  diameter  or  other  characteristic  dimension 
°F  =  temperature,  degrees  Fahrenheit 
Fd  =  drag  force 
Fg  =  force  of  gravity 
g  =  acceleration  of  gravity 
h  =  height  of  fall  at  time  t  or  velocity  v 
H  =  total  height  of  fall 

K  =  acceleration- coefficient  numeric  =  (p/g)  (A/m)  (u/D)2 


m,  M  =  mass 

m0  =  reference  mass 

P0  =  atmospheric  pressure,  ambient 

Rd  =  Reynolds  number  =  vD/y 

Rv  =  a  function  concerned  with  error  in  computed  velocity 
Ra  =  a  function  concerned  with  error  in  computed  a 
R'a  =  an  approximation  to  Ra 
R'^=  an  approximation  to  Ra 
s  -  frontal  area  of  the  object 
t  =  instantaneous  time 
t=  thickness  of  plate -like  objects 
t0  =  a  reference  thickness 
T  =  total  time  of  fall 

Tnd  =  total  time  of  fall  assuming  no  drag  forces 
AT  =  T  -  TNd 
6T  =  timing  error 
v  =  instantaneous  velocity 

V  =  impact  velocity 

V  nd  =  “no-drag”  impact  velocity 
AV  =  error  in  computed  impact  velocity 
X  =  time  numeric  s  tg(D/u) 

AXf=  increment  of  X 

Y  =  height-of-fall  numeric  =  hg(D/i/)2 
AY  =  increment  of  Y 

Z  =  velocity  numeric  =  Rd  =  vD/u 
a  =  acceleration  coefficient  =  sCj^/m 
cx0  =  reference  acceleration  coefficient 
Aa  =  error  in  computed  a 

eT  -  t/tnd 

9't  =  an  approximation  to  6? 

0%  =  an  approximation  to  6t 

eY  =  v/vND 
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0y  =  an  approximation  to  0V 
v  =  kinematic  viscosity  of  the  air 
p  -  density  of  air 
p  -  density  of  plate -like  objects 
p0  =  a  reference  density  of  plate -like  objects 
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Chapter  1 

INTRODUCTION 


1.1  OBJECTIVES 

During  the  1955  and  1957  Test  Operations  at  the  Nevada  Test  Site  (NTS),  masses  and 
velocities  of  more  than  20,000  objects,1-5  such  as  glass  fragments  from  windows,  stones,  steel 
fragments,  and  spheres,  which  were  energized  by  blast  winds  were  experimentally  determined, 
and  the  time -displacement  history  of  an  anthropometric  dummy  simulating  man5  was  measured. 
The  availability  of  such  a  mass  of  data  stimulated  an  analytical  study  calculated  to  arrive  at  a 
mathematical  formulation  capable  of  predicting  the  translation  of  objects  by  blast  winds.  The 
model7  designed  to  predict  this  translation  assumed  that  certain  aerodynamic  properties  of  the 
objects  displaced  were  known.  The  purpose  of  this  report  is  to  describe  the  drop-test  studies, 
both  theoretical  and  experimental,  which  were  undertaken  to  determine  the  aerodynamic 
properties  of  interest. 

1.2  HISTORICAL  BACKGROUND 

When  objects  dropped  in  air  were  first  observed  to  fall  at  varying  rates,  it  was  presumed 
that  the  difference  was  inherent  in  the  nature  of  the  materials  of  which  the  objects  were  com¬ 
posed.  In  the  early  16 th  century,  Galileo  disproved  this  theory  and  established  the  true  cause 
of  the  observed  variations.  Small  spheres  of  different  substances  were  dropped  at  the  same 
instant  from  the  Leaning  Tower  of  Pisa,  and  all  hit  the  ground  at  nearly  the  same  time.  Next, 
Galileo  dropped  objects  composed  of  the  same  materials  in  a  variety  of  shapes  and  observed 
that  they  then  fell  at  different  rates.  Galileo  concluded  that  gravity  acted  upon  all  objects  with 
the  same  intensity  but  that  the  air  resisted  their  fall  in  varying  degrees,  depending  on  the 
object  being  dropped.  It  then  followed  that  all  objects  would  fall  at  the  same  rate  in  a  vacuum. 

After  the  invention  of  the  air  pump,  Newton  verified  Galileo’s  theory  with  his  famous 
“guinea  and  feather”  experiment.  A  gold  coin,  called  a  “guinea,”  and  a  feather  were  placed  in  a 
tube  about  5  ft  long.  When  the  tube  was  turned  on  end,  they  were  observed  to  fall  at  different 
rates.  When  most  of  the  air  was  removed  from  the  tube  with  an  air  pump,  no  difference  in  the 
rates  of  fall  of  the  guinea  and  the  feather  was  perceivable. 

The  force  with  which  the  air  acts  upon  a  falling  object  came  to  be  known  as  a  drag  force 
since  it  always  resisted  the  motion  of  the  object.  The  first  experiments  to  determine  the 
nature  of  this  drag  force  were  performed  in  1853  by  Piobert,  Morin,  and  Didion,5  who  observed 
an  object  falling  in  air  and  recorded  its  velocity  until  it  reached  terminal  velocity.  Since  the 
object  at  that  point  was  not  undergoing  acceleration,  the  drag  force  equaled  the  force  of  gravity. 
They  were  able  to  obtain  terminal  velocities  up  to  about  30  ft/sec.  In  1892  Cailletet  and 
Colardeau9  dropped  plane  surfaces  of  various  shapes  from  the  Eiffel  Tower  and  recorded  the 
descent  times;  they  obtained  velocities  up  to  90  ft/sec.  In  1905,  Eiffel10  conducted  a  series  of 
experiments  that  involved  attaching  the  object  to  be  tested  to  an  apparatus  that  slid  down  a 


13 


stretched  wire.  This  apparatus  actually  recorded  both  the  distance  dropped  and  the  drag  force 
as  a  function  of  time,  measured  by  the  vibrations  of  a  tuning  fork.  Velocities  up  to  130  ft/sec 
were  obtained. 

Today  most  experiments  involving  drag  measurements  are  conducted  in  an  artificial  air 
stream.11  This  method  possesses  many  advantages  over  the  drop  test.  In  addition  to  improving 
accuracy,  the  use  of  an  artificial  air  stream  makes  it  easier  to  vary  the  environment  of  the 
object  being  tested.  The  drop  technique  is,  however,  still  used  to  determine  air  density,12-15 
to  study  parachute  descents,16-18  and  to  determine  the  drag  force  on  bombs  at  high  subsonic 
speeds.19  It  has  also  been  used  to  study  impact  damage  to  small  laboratory  animals.20-23 
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Chapter  2 


EQUATIONS  OF  MOTION 


2.1  INTRODUCTORY  REMARKS 

To  predict  the  motion  of  an  object  displaced  by  the  movement  of  the  surrounding  air,  one 
must  know  the  object’s  acceleration  coefficient  (a).  The  quantity  a  has  been  defined  as  the 
product  of  the  presenting  area  and  the  drag  coefficient  divided  by  the  mass  of  the  object.1  The 
same  quantity  is,  of  course,  involved  in  the  equations  of  motion  of  an  object  in  free  fall,  where 
the  object  is  moving  instead  of  the  air.  It  will  be  shown  in  this  chapter  by  development  of 
equations  of  motion  that  a  can  be  determined  experimentally  by  drop  tests  if  the  time  of  fall, 
height  of  drop,  acceleration  of  gravity,  and  air  density  are  known.  Analysis  of  experimental 
data  involves  the  assumptions  that  (1)  no  vertical  air  motion  existed,  (2)  the  air  density  re¬ 
mained  constant  over  the  entire  drop  interval,  and  (3)  the  a  of  the  falling  object  did  not  change 
during  the  fall.  Measures  were  taken  to  reduce  air  motion  during  the  drop  tests  (see  Chap.  3). 
Since  the  total  drop  height  was  less  than  50  ft,  the  assumption  of  constant  air  density  is  reason¬ 
able.  Changes  in  a  may  occur  during  fall  by  two  processes:  First,  the  presenting  area  may 
change  for  objects  other  than  spheres  because  of  rotation,  and,  second,  the  drag  coefficient,  a 
function  of  Reynolds  number  and  thus  of  velocity  of  fall,  may  change.  A  further  discussion  of 
these  topics  will  be  found  in  Chap.  4. 


2.2  BASIC  EQUATIONS 


A  free-falling  object  is  subjected  to  two  forces:  a  drag  force  (F^)  resulting  from  the  drag 
of  the  air  on  the  object  and  the  force  of  gravity  (Fg).  For  the  case  of  free  fall  in  still  air, 


Fd  =  (V2)PV2CDS 


(2.1) 


and 


Fg  =  mg 


where  Fd 

P 

v 

CD 

S 


Fg 

m 


g 


drag  force 
air  density 
velocity  of  object 
drag  coefficient 
frontal  area  of  object 
force  of  gravity 
mass  of  object 
acceleration  of  gravity 


(2.2) 
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These  two  forces  (Fd  and  Fg)  act  in  opposition.  By  Newton’s  second  law  of  motion,  their 
algebraic  sum  is  equal  to  the  product  of  the  mass  and  the  acceleration  (a)  of  the  object,  or 

ma  =  mg-(V2)pv2CDs 
which  can  be  written  in  the  form: 

<2-3’ 

The  quantity  sCD/m  is  the  acceleration  coefficient  (a).  Its  units  are  area/mass  (ft2 /lb  is 
used  throughout  this  report).  Thus 


Alpha  has  a  special  significance  since  it  combines  the  significant  parameters  of  the  object 
which  influence  the  drag  force  and  consequently  the  acceleration  of  the  object,  i.e.,  if  two 
objects  with  the  same  a  were  dropped  at  the  same  time  under  identical  conditions,  their  accel¬ 
erations  and  velocities  would  remain  the  same  throughout  the  fall. 

Therefore  Eq.  2.3  can  now  be  written  as 


By  definition,  v  and  a  can  be  written  in  differential  form: 
dh 


where  t  =  time  and  h  -  height. 
By  combining  Eqs.  2.6  and  2.7, 
v  dv 


and  substituting  a  from  Eq.  2.5  into  Eq.  2.8, 
v  dv 


1  2  pa\ 

gv 

‘-v  2i0 

Now,  if  a  is  assumed  to  be  constant  throughout  the  low- velocity  range  v*  <  v  <  v2,  Eq.  2.9 
can  be  integrated  (see  Appendix,  Sec.  A.  1.1)  over  the  range  hi  <  h  <  h2  and  tj  <  t  <  t2  to  give  the 
following  equations: 

h2  -  h4  =  (pa)-1  In  {[1  -  vt  (pa/ 2g)V2]  [1  +  vt  (pa/ 2g)Vl!]} 

-  ( pa)- 1  In  {[1  -  v2  (pa/2g)^]  [1  +  v2  (pa/2g)^]}  (2.10) 
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and 


t2  -  tj  =  (2g par*  In  {[1  -  v,  (pa/2g)^]  [1  +  v2  (pa/ 2g)v=]} 

-  (2gpa)-^2  In  {[1  -  v2  (pa/2g)‘/2]  [1  +  vt  (pa/2g)‘/2]}  (2.11) 

Before  these  equations  were  used  for  numerical  computations,  they  were  transformed  to  a 
more  convenient  nondimens ional  form,  which  is  explained  in  the  Appendix  (Sec.  A.  1.2).  Fig¬ 
ure  2.1  illustrates  how  the  drag  coefficient  of  a  sphere  varies  with  Reynolds  number.2  The 
Reynolds  number  (Rd)  is  a  dimensionless  quantity  defined  by 

vD 

Rd=—  (2.12) 

where  Rd  =  Reynolds  number 

v  =  velocity  of  the  sphere 
D  -  diameter  of  the  sphere 
v  =  kinematic  viscosity  of  the  air 

Figure  2.2  shows  the  Reynolds  number  as  a  function  of  velocity  and  characteristic  dimen¬ 
sion  (diameter  for  a  sphere,  length  of  side  for  a  square  plate,  etc.).  The  characteristic  dimen¬ 
sion  is  read  directly  in  inches;  whereas,  velocity  is  read  directly  only  if  the  kinematic  viscosity 
of  the  air  is  0.000179  ft2/sec.  This  value  was  chosen  because  it  represents  the  average  kine¬ 
matic  viscosity  at  5300  ft  above  sea  level,3  the  approximate  altitude  of  the  southeast  section  of 
Albuquerque,  N.  Mex.,  where  the  experiments  were  performed.  If  the  kinematic  viscosity  has 
any  other  value,  the  velocity  variable  must  be  modified  by  the  appropriate  factor,  as  indicated 
in  Fig.  2.2,  to  convert  it  to  its  true  velocity.  If  the  velocity,  diameter  of  a  sphere,  and  the  kine¬ 
matic  viscosity  of  the  surrounding  medium  are  known,  Fig.  2.2  can  be  used  to  find  the  Reynolds 
number.  Then  for  that  Reynolds  number,  the  corresponding  drag  coefficient  can  be  determined 
with  the  aid  of  Fig.  2.1. 


2.3  EQUATIONS  FOR  CONSTANT  a 

If  a  is  assumed  to  be  constant  throughout  a  fall;  if  the  initial  conditions  are  vA  =  ht  =  tj  =  0; 
and  if  v2  =  V,  h2  =  H,  and  t2  =  T  (see  Appendix,  Sec.  A.  1.1),  Eqs.  A. 3  and  A. 5  become 

V  =  (2gH)'/2  [(1  -  e~Pan)/p aH]H  (2.13) 

and 

T  =  [(2H/g)1/2 (pa H)  y2/2]  +  (2H/g)‘/2  (paH)~1/2  In  [1  +  (1  -  e_paH)1/2]  (2.14) 

An  alternative  way  of  writing  Eq.  2.14  is 

T  =  (2H/g)y2  (paH)-V2  cosh-1  e  P“H/2  (2.15) 


An  explicit  solution  of  Eq.  2.13  for  a  in  terms  of  V  and  H  was  not  possible;  however,  for  a 
fixed  value  of  H,  V  can  be  plotted  as  a  function  of  a.  Such  a  plot,  illustrated  in  Fig.  2.3,  was 
found  to  be  useful  in  a  previous  study  of  the  impact  loading  of  small  anesthetized  laboratory 
animals4  in  which  impact  velocity  was  determined  by  means  of  high-speed  photography.  The 
experimental  data  shown  on  this  chart  (Fig.  2.3),  however,  were  determined  in  the  present 
study  (Sec.  4.5). 

Acceleration  as  a  function  of  height  of  fall  can  be  found  by  substituting  V  from  Eq.  2.13 
into  Eq.  2.5  to  give 


a 


=  ge 


-paH 


(2.16) 
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2.4  EQUATIONS  WITHOUT  DRAG  * 


The  equations  for  free  fall  with  no-drag  force  and  initial  conditions  V  =  H  =  0  when  T  =  0 

are 

aND  =  g  (2-17) 

VND  =STND  (2.18) 

and 

H  =  g(TND)2/2  (2.19) 

Equation  2.19  can  be  solved  for  T  to  give 

Tnd  =  (2H/g)1/2  (2.20) 

Substituting  T  from  Eq.  2.20  into  Eq.  2.18  gives 

VND  =  (2gH)ya  (2.21) 

2.5  COMPARISON  OF  EQUATIONS  OF  MOTION  WITH  AND  WITHOUT  DRAG 
A  summary  of  the  preceding  equations  is  given  here: 


For  any  initial  conditions: 

(a)  With  drag  (for  constant  a) 

a=g(i-v2g) 

(b)  No  drag 

aND  =  g 

For  initial  conditions  v  =  h  =  0  when  t  =  0: 

(a)  With  drag  (for  constant  a) 

V  =  (2gH)1/2  [(1  -  e_pfflH)/paH]y2 
a  =  g  e~paH 

T  =  [(2H/g)1/2  (p«H)V2]  +  (2H/g)y2  (paH)~y2  In  [1  +  (1  -  e'^V2] 

(b)  No  drag 

VND  =  (2gH)y2 

aND  =  S 
Tnd  =  (2H/g)y2 


(2.5) 

(2.17) 


(2.13) 
(2.16) 

(2.14) 

(2.21) 

(2.17) 

(2.20) 


*The  subscript  ND  is  used  to  indicate  the  no-drag  value  of  the  parameter;  i.e.,  Tnd  is 
the  time  an  object  would  take  to  fall,  assuming  the  object  did  not  experience  drag  force. 
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For  the  initial  conditions  v  =  h  =  0  when  t  =  0,  the  quantities  with  drag  are  arranged  so 
that  they  differ  from  those  without  drag  by  a  factor  that  is  a  function  of  the  dimensionless 
quantity  pa H. 

Figure  2.4  shows  these  results  in  a  different  form.  Theta  (0)  is  defined  as  the  ratio  of  a 
variable  (V,  T,  or  a)- with  drag  to  the  same  variable  without  drag  and  is  a  dimensionless 
quantity  expressed  as  a  function  of  the  paH.  For  example,  0y  is  the  final  or  impact  velocity 
divided  by  the  velocity  for  a  drop  with  no  drag.  The  final  or  impact  velocity  can  be  computed 
from  Fig.  2.4  and  measured  values  of  H,  g,  and  T  (assuming  a,  g,  and  p  are  constant).  This 
calculation  can  be  made  by  computing  T^d  anc*  #T  an<3  then  using  the  0T  curve  to  determine 
paH.  For  that  value  of  paH,  is  then  read  and  multiplied  by  the  computed  value  of  Vnd  to 
determine  the  impact  velocity.  If  the  values  for  p  and  a  are  known  and  H  is  measured,  0 t  and 
0y  can  be  read  directly  from  Fig.  2.4. 

The  dotted  curves  in  Fig.  2.4,  which  are  approximations  to  the  solid  curves,  are  good  over 

certain  ranges  of  paH  and  are  useful  to  simplify  computations.  They  were  computed  from  the 

following  equations: 

9't  =  [(paH)H/2]  +  (paH)",/4  In  2  (2.22) 

and 

=  (paH)-1'4  (2.23) 

These  approximate  relations  depend  upon  the  fact  that  c~paU  —  0  very  fast  as  paH  —  °°.  Hence 
zero  was  substituted  for  in  the  equations  for  0x  and  0y  (Fig.  2.4)  to  arrive  at  the  ap¬ 

proximations.  Using  Eq.  A.16  from  the  Appendix  (Sec.  A. 2)  and  the  definition  0 v  =  V/Vnd  ,  one 
can  write 


V/Vt  =  9V  (2gH)1/a  (2g/pa)_^2 

=  0V  (paH)^  =  9V/9{,-  (2.24) 

where  Vtis  the  terminal  or  limiting  velocity  as  pa H  — -  Thus  the  convergence  of  0y  to  0y 
indicates  V  — -  Vt  or  the  approach  of  the  velocity  of  an  object  to  terminal  velocity  as  the 
height  of  fall  increases.  Another  approximation  to  9^, 

0"  =  1  +  (paH/12)  (2.25) 

for  low  values  of  paH  was  found  by  expanding  6> p  in  a  Maclaurin’s  series  as  follows: 

0T(paH)  =  6>t  lo  +  d@T/d(paH)  |0  (paH)  +  d2eT/d(paH)2  |0  (paH)2/2!+  . . . 

=  1  +  (paH/12)  +  (pa H)2/480  +  . . .  (2.26) 

The  approximation  shown  by  Eq.  2.25  is  then  the  first  two  terms  of  the  Maclaurin’s  series  for 
0T(paH). 

Although  impact  velocity  can  be  computed  from  height  of  fall  with  Fig.  2.4,  as  described 
in  the  preceding  paragraphs,  Fig.  2.5  represents  a  simpler  way  to  accomplish  the  same  thing. 
It  is  especially  useful  if  g  =  32.174  ft/sec2,  in  which  case,  V,  T,  and  H  can  be  read  directly. 
This  value  for  g  was  chosen  because  it  represents  the  acceleration  of  gravity5  at  sea  level.  If 
the  acceleration  of  gravity  has  any  other  value,  the  velocity  and  time  variables  must  be  modi¬ 
fied  by  the  appropriate  factors,  as  indicated  in  Fig.  2.5,  to  convert  them  to  actual  velocity  and 
time. 

2.6  ESTIMATIONS  OF  ERRORS 

To  estimate  the  percentage  of  error  in  computed  impact  velocity  (V)  introduced  by  a 
timing  error  (6T),  the  following  equation  was  used: 


20 


(2.27) 


100AV/V  =  100(AV / AT)(1/V)  AT  =  100(a/V)  6T 

Substituting  for  a  and  V  from  Eqs.  2.16,  2.13,  and  2.20  in  Eq.  2.27  gives 

100AV/V  =  [(paH)'/2e-P“H  (1  -  e-PaHf*]  1005T/TND 

=  RV  •  1005T/TND  (2.28) 


where  Ry  is  the  dimensionless  quantity  in  brackets.  Figure  2.6  is  a  plot  of  Rv  vs.  paH.  To 
find  the  estimated  percentage  of  error  in  computed  impact  velocity  due  to  error  in  measured 
time  of  fall,  determine  Rv  from  Fig.  2.6  and  use  this  value  in  Eq.  2.28.  It  is  interesting  to 
note  that  100AV/V  <  1006T/Tnd  for  all  values  of  paH;  indeed,  by  substituting 

TND  =  T{[(paH)^/2]  +  (paH)-y2  In  [l  +  (1  -  e-P“H)1/T1} 


from  Eqs.  2.14  and  2.20  into  Eq.  2.28,  it  can  be  shown  that  AV/V  <  ST/T  for  all  values  of  pa H. 
Thus  impact  velocity  can  be  determined  with  greater  accuracy,  percentagewise,  than  the 
measurement  of  time  of  drop  from  which  it  is  computed. 

A  similar  procedure  was  used  to  estimate  the  percentage  of  error  in  computed  a  intro¬ 
duced  by  a  timing  error  6T. 


100A  a/a  =  100(Aa/AT)(l/a)AT  =  100 


<5T 


(2.29) 


Differentiate  Eq.  2.14  with  respect  to  a  and  substitute  in  Eq.  2.29  to  obtain 

100A a/a  =  2{[(paH)VV2]  -  (paH)”*  In  [1  +  (1  -  e^11)*]  +  (paU)'A  e-P«H[(i  _  e-pa H) 

+  (1  -  e-f’aH)1''*]-1}-1  1006T/TND 

=  Ra  1006T/Tnd  (2.30) 

where  the  dimensionless  quantity  in  braces  is  given  the  symbol  Ra. 

Figure  2.7  is  a  plot  of  Ra  vs.  pa H.  To  find  the  estimated  percentage  of  error  in  com¬ 
puted  a  resulting  from  error  in  measured  time  of  fall,  determine  Ra  from  Fig.  2.7  and  use 
this  value  in  Eq.  2.30.  The  two  dashed  curves  are  reasonable  approximations  to  the  solid 
curve  over  certain  ranges  of  pa H.  The  R^  approximation  was  formed  in  the  same  way  the  6j 
and  approximations  were  formed  (Eqs.  2.22  and  2.23).  The  R^  approximation  for  low 
values  of  pa H  was  obtained  by  differentiating  Eq.  2.25  and  substituting  it  into  Eq.  2.29  to  obtain 

R'a  =  12/paH  (2*31) 


Equation  2.31  was  used  to  approximate  the  timing  accuracy  needed  for  the  experiment, 
where  p  «  0.06  lb/ft3,  H  «  50  ft,  g  «  32  ft/sec2,  and  a  «  0.03  ft2/ lb.  Since  the  timing  require¬ 
ment  is  greatest  for  low  values  of  a,  the  value  used  in  this  report  was  the  lowest  one  antici¬ 
pated  for  the  objects  to  be  dropped.  From  Eq.  2.31,  Ra  «  12/paH  =  12/0.09.  If  this  value  of 
Rq,  is  used  for  the  quantity  in  braces  in  Eq.  2.30, 

100Aa/a  =  (12/0.09)1006T[2(50)/32]“^  «  75006T 

Thus,  if  6T  =  0.2  msec,  then  lOOAa/a  «  1.5  per  cent.  Therefore  to  limit  the  error  in  com¬ 
puted  a  to  about  1  per  cent,  the  timing  mechanism  would  need  to  determine  the  drop  time  to 
the  nearest  0.1  msec. 
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Fig.  2.3  —  Impact  velocity  vs.  height  of  fall  and  a.  Included  in  the  figure  are  the  animal  data  from  the 
present  report  showing  mean  values  (represented  by  o)  and  standard  deviations  (represented  by  *  ’ ) 
of  the  experimentally  determined  a . 
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Find  Rv  from  curve  and  percentage  of  error  in  V  from  equation  shown. 


Fig.  2.7 — Method  for  estimating  error  in  computed  a  resulting  from  error  in  measured  time  of  fall 
Find  Rr  from  curve  and  percentage  of  error  in  a  from  equation  shown. 


Chapter  3 


EXPERIMENTAL  PROCEDURE 


The  experimental  system  used  to  measure  the  time  required  for  various  objects  to  fall  a 
measured  distance  is  shown  diagrammatically  in  Fig.  3.1.  The  object  to  be  dropped  was  placed 
near  the  end  of  a  bar  that  was  designed  to  swing  about  its  mid -point  (see  drop  device  in  upper 
right  of  Fig.  3.1).  The  drops  were  initiated  manually  by  tripping  the  catch  shown  just  left  of 
the  microswitch.  This  action  allowed  a  helical  spring  to  rotate  the  bar  about  its  mid-point.  As 
the  bar  began  to  move,  the  test  object  was  left  free  to  fall;  and,  at  the  same  time,  the  micro¬ 
switch  was  activated,  starting  an  electronic  timer  (Hewlett  Packard  Electronic  Counter,  model 
522 -B).  At  the  termination  of  the  fall,  the  test  object  struck  an  aluminum  plate,  creating  a 
noise,  which  was  detected  by  a  contact  microphone  placed  on  the  underside  of  the  impact  plate. 
The  signal  from  the  microphone  then  stopped  the  electronic  timer. 

Figure  3.2  is  a  picture  of  the  drop  apparatus  and  a  close-up  view  of  the  release  mecha¬ 
nism.  The  dropping  bar,  supported  on  both  sides  by  vertical  members,  was  made  thin  (3/s  in.) 
to  minimize  air  disturbance  during  release  of  the  test  object.  Plates  of  various  sizes  were 
made  to  be  attached  to  the  right  end  of  the  dropping  bar,  whose  total  length  was  32  in.,  so  that 
objects  as  large  as  a  rabbit  could  be  conveniently  dropped.  A  helical  spring,  hidden  from  view 
by  the  vertical  support,  was  designed  to  ensure  that  the  downward  acceleration  of  the  right 
side  of  the  bar  would  be  greater  than  the  acceleration  of  gravity.  After  the  bar  moved  approxi¬ 
mately  90°,  its  energy  was  partially  dissipated  by  impact  with  a  small  pad  of  Ensolite  (a 
product  of  the  U.  S.  Rubber  Co.)  placed  between  the  upright  supports;  thus  rebound  was  less¬ 
ened,  preventing  interference  of  the  bar  with  the  falling  object. 

It  was  necessary  to  pay  particular  attention  to  the  design  and  construction  of  the  release 
mechanism.  If  any  motion  had  been  imparted  to  the  dropping  bar  during  the  process  of  release, 
then  the  times  of  drop  would  not  have  been  reproducible.  Therefore,  the  sliding  parts  of  the 
latch  were  made  of  blocks  of  tool  steel,  precision  ground  and  clamped  at  the  ends  of  the  drop¬ 
ping  bar;  and  the  vertical  bar,  shown  on  the  left  side  of  Fig.  3.2,  was  designed  to  pivot  about  a 
point  near  the  base  of  the  apparatus.  Even  after  taking  these  precautions  in  design,  it  was 
found  that  more  reproducible  times  of  drop  could  be  obtained  if  the  catch  were  released  slowly 
by  means  of  the  knurled  nut  shown  at  the  left  in  Fig.  3.2,  thus  minimizing  motion  induced  in 
the  dropping  bar. 

Tests  were  conducted  in  an  elevator  shaft  in  which  objects  were  dropped  48.372  ft.  Pre¬ 
cautions  were  taken  to  minimize  vertical  air  currents  by  keeping  all  access  doors  to  the  shaft 
closed. 

The  impact  plate  shown  schematically  in  Fig.  3.1  consisted  of  a  V4-in. -thick  2-ft2  alumi¬ 
num  plate.  The  plate  was  fastened  to  a  table-like  support  about  6  in.  high  made  of  3/4-in. -thick 
plywood  in  which  a  hole  was  made  in  the  center  to  accommodate  the  contact  microphone.  It  is 
pertinent  to  note  that  the  impact  plate  made  a  relatively  small  target  for  irregular  objects 
dropping  48  ft.  Occasional  “misses”‘were  experienced,  especially  by  the  smaller  asymmetri¬ 
cal  objects.  Plate -like  objects  dropped  flat  also  tended  to  veer  away  from  the  impact  plate  if 
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a  small  degree  of  rotation  occurred.  Surprisingly  enough,  the  release  operator  could  usually 
observe  the  object’s  fall  well  enough  to  detect  the  presence  or  absence  of  rotation.  Thus  the 
data  included  in  this  report  are  for  objects  that  had  little  or  no  rotation  while  falling. 


MICROPHONE 


Fig.  3.1 _ Block  diagram  of  the  drop-test  equipment  showing  release  apparatus  and  starting 

and  stopping  circuits. 
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Fig.  3.2  —  Release  apparatus  used  in  the  drop  tests.  The  helical  spring  designed  to  accelerate  the  motion  of  the  dropping 
bar  (shown  in  horizontal  position)  is  hidden  from  view  behind  the  vertical  support.  The  inset  is  an  enlargement  of  the 
release  mechanism  (see  text). 


Chapter  4 


ANALYSIS  OF  DATA 


4.1  CONSTANTS  FOR  THE  DROP  TESTS 
The  constants  for  the  drop  tests  were 
g  =  32.126  ft/sec2 
H  =48.372  ft 


p  =  0.060437  lb/ft3 
v  =  0.0002071  ft2  /sec 

The  value  of  g  was  obtained  from  the  Coast  and  Geodetic  Survey,  Department  of  Com¬ 
merce.1  H  was  measured  directly;  whereas,  p  and  v  were  computed  by  measuring  the  air 
pressure  and  temperature  and  substituting  them  in  the  following  derived  equations:2 


2.6954  Pq  (psi) 
459.69  +  t(°F) 


2.7803  (10)-7  [456.69  +t(°F)] 


(4.1) 

(4.2) 


Since  the  drop  tests  were  conducted  over  a  period  of  three  days,  the  above  values  of  p 
and  v  are  averages;  however,  the  deviations  about  these  average  values  from  day  to  day  were 
small,  the  range  in  both  cases  being  less  than  1  per  cent. 


4.2  TIME  CORRECTIONS 

Steel  spheres,  7/i6  and  Vi6  in-  in  diameter,  were  dropped  during  each  of  the  three  days  when 
drop  tests  were  conducted.  Although  the  drop  times  for  a  particular  sphere  on  a  given  day 
were  quite  consistent  (low  standard  error),  the  mean  drop  time  varied  significantly  from  day 
to  day.  Since  the  air  density  remained  nearly  constant,  the  differences  were  probably  due  to 
minor  adjustments  that  had  been  made  on  the  microswitch  and  the  release  apparatus  between 
daily  tests. 

The  circumstance  stated  above  pointed  to  the  need  for  accurately  computed  drop  times 
for  the  two  “calibration”  spheres  in  order  to  determine  the  timing  deviations  inherent  in  the 
experimental  equipment.  This  was  accomplished  by  the  stepwise  integration  of  Eqs.  2.10  and 
2.11,  the  height  of  drop  and  time  of  drop  being  computed  for  small  steps  in  velocity  for  which 


CD  (and  thus  a)  could  be  considered  to  have  a  constant  value  (see  Fig.  2.1).  Thus  the  following 
values  were  computed  for  time  above  no -drag  time  (AT):  For  a  3/16 -in. -diameter  steel  sphere, 
AT  =  35.05  msec;  and  for  a  7/16 -in. -diameter  steel  sphere,  AT  =  14.91  msec. 

These  values,  when  compared  to  the  measured  times  of  drop,  indicated  that  the  timing 
error  for  each  day  was  essentially  the  same  for  the  two  spheres.  The  average  corrections  to 
be  added  to  the  measured  times  of  fall  were  calculated  to  be:  first  day,  +3.45  msec;  second 
day,  -0.14  msec;  and  third  day,  -1.65  msec. 

In  addition  to  these  corrections,  another  correction  was  needed  for  the  tests  when  heavy 
objects  were  dropped  on  a  pad  of  Ensolite.  This  1 -in. -thick  pad  was  used  on  the  third  day  to 
help  protect  the  contact  microphone.  Several  spheres  were  dropped  both  with  and  without  the 
pad.  The  differences  in  times  of  fall  produced  by  the  pad  were  plotted  against  the  uncorrected 
times  of  fall  with  the  pad.  The  points  indicated  a  linear  relation.  This  correction  was  then 
combined  with  the  constant  correction  previously  determined  for  the  third  day  to  give  a  total 
time  correction  (Fig.  4.1). 

4.3  DATA  FOR  SPHERES,  STONES,  AND  MILITARY  DEBRIS 

4.3.1  Introductory  Remarks 

It  was  implied  in  Sec.  4.2  that  to  determine  an  accurate  time  of  drop  relation  for  spheres, 
it  was  necessary  to  allow  for  a  varying  a  dependent  on  the  Reynolds  number,  which  changes 
throughout  the  fall.  If  this  were  true  for  irregular  objects,  then  the  drop-test  technique  would 
be  doomed  to  failure  since  the  relation  between  drag  coefficient  and  Reynolds  number  is  not 
known  for  objects  of  interest  to  this  study.  In  this  section  it  will  be  shown  that  the  dependence 
of  a  upon  Reynolds  number  for  irregular  stones  is  evidently  quite  different  from  that  for 
spheres.  This  difference  makes  it  reasonable  to  assume  for  experimental  purposes  that  a  for 
stones  can  be  considered  to  be  constant  for  the  entire  drop  interval. 

4.3.2  Sphere  Data 

First,  it  will  be  necessary  to  examine  the  test  data  for  the  spheres.  Figure  4.2  shows  a 
plot  of  “theoretical”  a  (based  on  CD  =  0.47)  for  spheres  vs.  measured  time  of  drop  in  excess 
of  the  no -drag  time.  The  curve  shown  was  computed  assuming  a  constant  a  by  using  Eq.  2.14. 
For  most  of  the  spheres,  the  times  measured  were  either  too  low  or  too  high  when  compared 
to  the  constant  a  curve.  An  examination  of  the  lower  curve  in  Fig.  4.3  will  help  explain  those 
deviations.  Marked  on  this  chart  are  the  Reynolds  numbers  computed  for  the  final  or  impact 
velocities  for  various  spheres.  The  position  of  the  Vig- in. -diameter  steel  sphere,  for  instance, 
indicates  that  high  drag  coefficients  existed  for  a  large  percentage  of  the  fall.  The  effect  of 
this  was  that  the  measured  time  of  drop  was  longer  than  that  expected  assuming  a  constant 
drag  coefficient  of  0.47  (Fig.  4.2).  Times  for  the  % -in. -diameter  steel  spheres,  on  the  other 
hand,  corresponded  reasonably  well  with  the  constant  a  assumption.  Undoubtedly  the  reason 
for  this  was  that  the  effect  of  the  initially  high  drag  coefficients  were  compensated  for  by  co¬ 
efficients  lower  than  0.47  near  the  termination  of  the  drop3  (Fig.  4.3).  For  spheres  still  larger 
in  diameter  (e.g.,  3/4-in. -diameter  aluminum  spheres),  the  drop  times  were  less  than  expected 
since  a  larger  percentage  of  the  drop-interval  involved  low  drag  coefficients  (Fig.  4.3).  The 
sphere  drop  involving  the  highest  Reynolds  numbers  was  that  of  a  3.34-in. -diameter  croquet 
ball  whose  weight  had  been  increased  with  a  brass  plug.  In  this  instance  a  large  portion  of  the 
drop  entailed  drag  coefficients  whose  values  were  near  0.47.  Thus  the  measured  drop  time 
corresponded  closely  to  that  expected  for  a  constant  coefficient  of  this  value.  It  is  interesting 
to  note  that  the  “effective”  drag  coefficient  for  the  entire  drop  interval  can  be  approximated  in 
each  individual  instance  by  the  coefficient  corresponding  to  one -half  the  Reynolds  number. 

4.3.3  Stones  Compared  to  Spheres 

To  study  the  differences  in  the  Reynolds  number -drag  relation  between  spheres  and 
irregular  objects,  it  was  necessary  to  analyze  the  two  groups  of  data  similarly.  This  involved 
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computing  effective  a  for  each  sphere  drop,  assuming  that  the  drag  coefficient  was  constant  for 
the  entire  drop  (Eq.  2.14).  The  a  values  so  computed  were  then  plotted  as  a  function  of  sphere 
mass,  as  illustrated  in  Fig.  4.4.  If  the  drag  coefficient  is  assumed  to  be  constant,  it  can  be 
shown  that  a  for  spheres  of  a  constant  density  varies  inversely  as  mass  to  the  V3  power.  Thus, 
if  the  drag  coefficient  had  the  same  constant  value  for  all  sizes  of  like-density  spheres,  the 
least-squares  regression  lines  shown  in  Fig.  4.4  would  have  slopes  of  —  /$.  Deviations  from 
this  value  indicate  that  the  smaller  spheres  had  effective  drag  coefficients  somewhat  greater 
than  those  for  the  larger  spheres. 

Statistical  data  pertinent  to  the  slopes  of  the  least-squares  regression  lines  for  spheres 
and  stones  are  recorded  in  Table  4.1.  In  addition  to  the  slopes,  the  standard  deviations  of  the 
slopes  were  computed;  and  from  these  data  the  probability  (pr)  that  the  measured  slopes  could 
have  been  deviations  from  the  “true”  slope  of  —  */3  was  computed.  Although  the  probabilities  for 
the  individual  cases  varied  considerably,  that  for  a  combination  of  all  spheres  proved  to  be  only 
0.034.  Stated  in  physical  terms,  the  Reynolds  number-drag  relation  for  spheres  is  such  that, 
for  the  purposes  of  the  present  study,  the  effective  drag  coefficient  varies  from  one  size  sphere 
to  another. 

Table  4.1  also  presents  the  results  of  a  similar  statistical  treatment  for  six  types  of  ir¬ 
regular  stones.  With  one  exception,  pumice,  the  probabilities  that  the  measured  slopes  could 
have  been  deviations  from  the  “true”  slope  of  -%  were  quite  high.  This  finding  indicated  that 
the  computed  slopes  are  likely  to  be  random  variations  from  the  true  slope.  The  variation  of 
the  pumice  data  from  the  -%  slope  was  probably  due  to  the  absorption  of  moisture  from  the 
air  in  the  outer  layer.  If  this  layer  were  of  uniform  thickness,  the  smaller  stones  would  suffer 
a  larger  relative  increase  in  weight  than  would  the  larger  stones. 

4.3.4  Stone  Data 

The  reasoning  presented  above  seemed  to  justify  a  constant-a  assumption.  The  a  values 
for  individual  objects  were,  therefore,  determined  graphically  from  a  curve  of  a  vs.  AT  similar 
to  that  shown  in  Fig.  4.2. 

Figure  4.5  is  a  plot  of  experimentally  determined  a  for  41  pieces  of  painted  gravel  as  a 
function  of  mass.  These  stones  were  taken  from  a  sample  of  stones  that  was  marked  for  identi¬ 
fication  and  used  in  secondary-missile  experiments  in  Nevada  during  the  Operation  Plumbbob 
weapons  tests.4  Since  the  stones  actually  used  in  these  experiments  had  masses  less  than  1500 
mg,  a  separate  analysis  was  made  for  these.  In  addition,  another  analysis  was  made  for  all 
stones,  assuming  a  slope  of  the  least-squares  regression  line  of  -V3  (Sec.  4.3.3).  The  results 
of  both  analyses  are  shown  in  Fig.  4.5. 

Drop  tests  were  made  with  limestone  fragments  so  that  there  could  be  a  comparison  of  the 
results  with  stones  whose  densities  were  less  uniform.  The  slope  of  the  least-squares  regres¬ 
sion  line  determined  for  these  data  (Fig.  4.6)  is  reasonably  near  the  expected  value  of  -% 

The  data  shown  in  Figs.  4.7,  4.8,  and  4.9  are  for  stones  obtained  at  three  different  loca¬ 
tions  (10P4,  1S1A,  and  4.3  GTS)  at  NTS  during  Operation  Plumbbob.4  An  analysis  was  made  of 
the  combined  data  for  stones  at  the  10P4  and  1S1A  locations  because  the  individual  analyses 
were  remarkably  similar. 

Figure  4.10  contains  the  data  for  the  pumice  fragments  that  were  discussed  in  Sec.  4.3.3. 
The  least-squares  regression  lines  for  all  the  stone  data  are  displayed  in  Fig.  4.11.  An  exami¬ 
nation  of  this  figure  points  to  the  similarity  between  pumice  and  stones  at  the  4.3  GTS  location, 
both  having  higher  a  for  given  stone  masses  than  the  other  samples. 

4.3.5  Military-debris  Data 

The  military  debris  used  in  the  drop  tests  consisted  of  steel  shrapnel  and  a  few  other  ir¬ 
regular  pieces  of  steel,  all  of  which  had  been  used  and  recovered  as  secondary  missiles  in  the 
Operation  Plumbbob  study.4  The  a-mass  data  for  these  steel  fragments  are  shown  in  Fig. 

4.12.  For  comparison  purposes,  the  theoretical  o'-mass  relation  for  steel  spheres,  assuming 
CD  =  0.47,  is  shown  on  the  same  chart.  It  is  apparent  and  reasonable  that,  for  a  given  mass,  a 
fragment  of  steel  has  a  much  larger  a  than  does  a  steel  sphere.  The  slope  of  the  least-squares 
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regression  line  is  less  than  that  for  the  theoretical  line  for  spheres;  however,  it  is  not  as  flat 
as  a  line  for  plates  (Sec.  4.4). 


4.4  DATA  FOR  PLATES,  CUBES,  AND  GLASS  FRAGMENTS 

4.4.1  Introductory  Remarks 

Drag  properties  of  window-  and  plate -glass  fragments  were  of  interest  to  this  study  be¬ 
cause  they  constitute  biological  hazards  when  hurled  through  the  air  by  the  action  of  a  blast 
wave  on  a  structure  containing  windows.4  Regular  plates  and  cubes  were  studied  only  for  the 
purpose  of  better  interpretation  of  the  glass -fragment  data  obtained  in  full-scale  nuclear  tests. 

4.4.2  Regular  Plates  and  Cubes 

It  should  be  pointed  out  that  the  plates  referred  to  in  this  study  are  really  “thick”  plates, 
i.e.,  no  attempt  was  made  to  sharpen  the  edges  in  order  to  achieve  the  “thin”  plate  that  is 
usually  used  in  wind-tunnel  studies.  For  this  reason  the  Reynolds  number- drag  relation 
shown  in  the  upper  part  of  Fig.  4.3  does  not  strictly  apply  to  the  experimental  data  reported 
here.  However,  it  is  useful  to  note  in  Fig.  4.3  that  within  the  region  of  interest  the  variation  in 
drag  coefficient  is  less  than  5  per  cent  (from  1.13  at  Rd  =  2  x  104  to  1.17  at  Rd  =  8  x  104). 

Figure  4.13  presents  drag  coefficients  for  regular  plates  in  normal  flow  determined  from 
drop-test  data,  assuming  that  a  was  constant  throughout  the  fall.  If  the  plates  used  had  been 
thin  plates,  those  with  the  smallest  area  could  have  been  expected  to  have  the  highest  drag 
coefficients  (Fig.  4.3).  Since  the  reverse  was  true,  it  can  be  concluded  that  the  thickness  effect 
must  outweigh  the  effect  of  Reynolds  number.  It  is  interesting  to  note  that  no  difference  in  CD 
could  be  detected  between  plates  that  were  square,  circular,  or  triangular  (equilateral). 

Similar  data  are  presented  in  Fig.  4.14  for  cubes  made  of  various  materials.  Although  a 
considerable  change  in  terminal  Reynolds  number  from  the  smallest  to  the  largest  cubes  was 
noticed,  no  significant  trend  in  the  drag  coefficient  data  was  apparent. 

With  information  obtained  from  the  same  drop  tests  for  cubes,  a  was  plotted  as  a  function 
of  the  length  of  the  cube  on  logarithmic  paper  (Fig.  4.15).  If  drag  coefficient  were  constant, 
then  the  following  relation  can  be  stated  for  cubes  made  of  material  of  constant  density: 
a  ~  (s/M)  ~  (length  of  side)^1.  Thus  a  comparison  of  the  data  points  in  Fig.  4.15  to  the  lines 
drawn  with  a  slope  of  - 1  indicate  the  degree  to  which  drag  coefficient  is  constant  for  the 
different  types  of  cubes. 

4.4.3  Data  for  Glass  Fragments 

Two  types  of  glass  fragments  were  used  in  this  study:  those  originating  from  window 
panes  (0.125  in.  thick)  and  those  from  plate-glass  panes  (0.225  in.  thick).  Figure  4.16  illus¬ 
trates  the  random  shapes  of  typical  window-glass  fragments.  The  number  appearing  with  each 
fragment  indicates  its  mass  in  grams. 

Drop  times  obtained  for  380  window-glass  fragments  were  used  to  determine  individual  a 
by  the  method  previously  described,  assuming  constant  a  throughout  the  fall.  When  a  was 
plotted  as  a  function  of  mass  on  logarithmic  paper,  it  was  apparent  that  the  points  tended  to 
form  not  one,  but  several  line  segments.  The  following  procedure  was  used  to  determine  the 
mass  boundaries  for  each  line  segment.  The  data  were  divided  into  small  mass  groups,  and 
the  mean  a  and  standard  deviation  in  a  were  determined  for  each  group.  These  data  were 
plotted  as  illustrated  in  Fig.  4.17.  The  mass  limits  for  each  line  segment  were  then  estimated, 
and  these  limits  used  to  determine  the  least-squares  regression  equations  for  each  segment. 
These  segments  are  shown  in  Fig.  4.18,  along  with  the  data  points  used  in  their  determination. 
Because  the  smaller  fragments  were  approximately  spherical,  the  regression  line  for  them 
was  given  a  theoretical  slope  of  —  */3  (Sec.  4.3.3).  Another  group  of  window  glass  was  dropped 
edgewise  (shown  as  triangles  in  Fig.  4.18).  The  regression  line  for  these  data  passed  near  the 
intersection  of  the  first  two  segments.  For  the  sake  of  simplicity,  another  regression  line 
was  computed  for  the  data,  with  the  requirement  that  it  pass  through  the  intersection  of  the 
first  two  segments.  The  standard  error  of  estimate  was  increased  very  little  by  this  procedure. 
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Alpha  data  for  plate  glass  (0.225  in.  thick)  dropped  flat  are  presented  in  Fig.  4.19.  It  is 
interesting  to  note  that  a  larger  scatter  is  shown  in  the  data  for  small  fragments  than  for  large. 
This,  no  doubt,  can  be  attributed  to  the  greater  mechanical  stability  during  fall  of  large  plates 
than  of  small  ones.  This  scatter  also  was  observed  at  the  time  of  the  experiments.  The  same 
distribution  of  scatter  can  be  detected  in  the  data  for  window  glass  dropped  flat  as  illustrated 
in  Fig.  4.18. 

A  summary  of  all  glass -fragment  data  is  shown  in  Fig.  4.20.  The  data  represented  by  line 
segment  1  for  the  smaller  fragments  can  be  used  for  either  plate  or  window  glass  in  any 
orientation  since  the  fragments  were  small  enough  that  their  flat  dimensions  were  not  ap¬ 
preciably  different  from  their  thickness  dimensions.  Thus  for  plate  glass  line  segment  1  was 
extended  to  higher  masses  until  it  intersected  line  segment  6  for  plate-glass  fragments  dropped 
flat.  This  extension  was  possible  because  the  densities  of  the  window  and  plate  glass  were  the 
same. 

4.4.4  Analysis  for  Thick  Plates  of  Random  Shape 

The  purpose  of  this  section  is  to  derive  relations  that  make  it  possible  to  determine  drag 
coefficient  (and  thus  a)  for  plates  of  any  thickness  and  presenting  area  within  the  limits  of  the 
experimental  data.  The  results  apply  to  plates  of  random  shapes  such  as  those  which  have 
been  described  (see  Fig.  4.16).  The  assumption  is  made,  and  verified  by  the  experimental 
data,  that  plates  with  the  same  characteristic  length -to -thickness  ratio  have  the  same  drag 
coefficient,  at  least  for  the  experimental  conditions  being  considered  here. 

The  characteristic  length  of  a  plate  of  random  shape  was  taken  to  be  the  length  of  a  side  of 
a  square  plate  of  an  equivalent  area,  i.e.,  the  square  root  of  the  area. 

It  has  been  shown  (see  Figs.  4.19  and  4.20)  for  wide  ranges  of  masses  and  for  plates  of 
constant  thickness  and  density  that  the  relation  between  the  acceleration  coefficient  and  mass 
can  be  expressed  by  an  equation  of  the  form 

log  a0  =  J  +  B  log  m0  (4-3) 

where  J  and  B  are  constants.  With  this  expression,  the  relation  between  Cq  and  s  Vt,  where  s 
represents  the  area  of  a  plate  in  normal  flow  and  t  represents  its  thickness,  can  be  determined. 

The  subscript  zero  is  used  to  designate  the  parameters  of  the  plates  used  in  this  study. 

The  same  terms  without  subscripts  represent  quantities  that  can  vary  between  limits  deter¬ 
mined  by  the  experimental  data.  Thus  a,  the  acceleration  coefficient,  may  be  defined  as  a0  = 
s0CD/m0  and  a  =  sCd/hi  (Sec.  2.2),  where  the  drag  coefficients  (Cd)  are  identical  for  the  con¬ 
dition  that  s01/2/to  =  sVt.  For  each  value  of  s  Vt  and,  therefore,  for  each  value  of  CD,  it  can  be 
stated  that  s  =  k4t2p,  s0  =  k^Jpo,  m  =  k2t3p,  and  m0  =  k2tjjp0,  where  k*  and  k2  are  constants,  t  is 
the  thickness  of  the  plate,  and  p  is  the  density  of  the  material  from  which  the  plate  is  made.  If 
these  relations  are  used,  the  following  equations  can  be  derived:  m0  =  mp0(t0/t)3  and  a0  = 
a(pt/p0t0).  By  the  use  of  these  equations  to  eliminate  m0  and  a0  from  Eq.  4.3,  the  following 
equation  is  obtained: 

log  [a(tp/t0Po)]  =  J  +  B  log  [m(t0/t)3  (p0/p)]  (4-4) 

Since  m  =  tps  and  a  =  sCp/m,  Eq.  4.4  can  be  reduced  to 

log  (CD/t„po)  =  J  +  B  log  [p0tjj(sVt)2]  (4.5) 

Since  t0  and  p0  were  assumed  to  have  constant  values  in  the  analysis  expressed  by  Eq.  4.3, 

Eq.  4.5  can  be  simplified  to 

log  CD  =  K,  +  K2  log  (s  Vt)  (4.6) 

where  Kt  and  K2  are  constants.  Thus  log  CD  is  a  linear  function  of  log  s  Vt  if  log  a  can  be 
assumed  to  be  a  linear  function  of  log  m  (Eq.  4.3). 
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To  verify  the  relation  derived  above,  a  plot  was  made  of  Cp  vs.  s1/2/t,  using  all  the  data 
for  plates  obtained  in  this  study.  Figure  4.21  shows  this  plot,  along  with  the  results  of  a  least- 
squares  regression  analysis.  The  data  presented  here  are  for  40  regular  plates,  referred  to 
in  Fig.  4.13,  and  77  irregular  plates,  referred  to  in  line  segments  4  and  6  in  Fig.  4.20.  It  is 
instructive  to  note  that  the  drag  coefficients  for  regular  plates  tend  to  be  slightly  higher  than 
those  for  the  irregular  ones.  Nevertheless,  the  standard  error  of  estimate  in  CD  for  all  plates 
was  only  approximately  5.6  per  cent,  being  between  that  for  window-glass  fragments  (3.4  per 
cent)  and  that  for  plate-glass  fragments  (7.8  per  cent).  The  areas  for  both  types  of  glass 
fragments  were  computed  from  mass,  thickness,  and  the  measured  glass  density  (39.657  g/in.3). 

4.5  SMALL -ANIMAL  DATA 

Six  mice,  six  rats,  six  guinea  pigs,  and  five  rabbits  were  dropped  during  these  tests.  The 
animals,  which  had  been  given  lethal  doses  of  Nembutal  before  the  drops,  were  dropped  as 
many  as  five  times  each. 

The  animals  were  dropped  in  two  different  initial  orientations:  first,  with  one  side  toward 
the  direction  of  fall  (lateral)  and,  second,  with  the  abdomen  toward  the  direction  of  fall  (ven¬ 
tral).  Figure  4.22  is  a  photograph  of  the  presented  areas  of  a  few  of  the  animals  in  these 
positions.  The  areas  of  most  of  the  animals  were  measured  from  similar  photographs  and  are 
recorded  in  Fig.  4.23,  with  circles  representing  the  lateral  areas  and  squares  representing 
the  ventral  areas.  The  straight  line  on  this  plot  has  a  slope  of  2/3.  It  tends  to  show  that,  at 
least  for  the  mice,  rats,  and  guinea  pigs,  the  area  is  proportional  to  the  mass  to  the  2/3  power 
(see  Ref.  5).  It  should  be  remembered  that,  although  the  animals  tended  to  hold  their  initial 
orientations  as  they  fell,  the  areas  on  Fig.  4.23  are  larger  than  those  “seen”  by  the  wind  since 
the  hair  tends  to  lie  close  to  the  body  as  the  animal  falls. 

Figure  4.24  is  a  plot  of  the  experimentally  determined  a  as  a  function  of  mass;  the  circles 
represent  the  a  of  animals  dropped  laterally  and  the  squares  represent  the  a  of  animals 
dropped  ventrally.  The  regression  line  was  computed  using  all  the  points.  The  slope  of  -0.3240 
is  near  the  theoretical  slope  of  -  which  was  determined  by  assuming  a  constant  density  and 
a  constant  drag  coefficient  for  any  group  of  similar  objects  (see  Secs.  4.3.3  and  4.3.4). 

Figure  2.3  shows  the  impact  velocity  as  a  function  of  drop  height  and  a.  The  a  values, 
determined  experimentally  from  the  animal  data,  are  plotted  on  the  theoretical  line  for  the 
actual  height  of  drop.  In  addition,  the  means  and  standard  deviations  of  the  experimentally  de¬ 
termined  a  are  tabulated  in  Fig.  4.24.  They  are  shown  in  Fig.  2.3  and  can  be  used  to  estimate 
the  means  and  standard  deviations  of  the  impact  velocities  of  the  various  animals. 

4.6  DRAG  PROPERTIES  FOR  MAN 

The  least-squares  regression  line  in  Fig.  4.24  can  be  extrapolated  up  to  a  mass  of  168  lb 
(75,977  g)  to  give  an  estimation  of  the  a  for  man  (0.0269  ft2/lb). 

The  extrapolated  value  for  the  a  for  man  was  compared  to  Schmitt’s  report  on  wind- 
tunnel  investigations  of  air  loads  on  human  beings6  whose  average  weight  was  168  lb.  The  data 
from  this  report  were  used  to  compute  and  average  the  a  for  a  group  of  eight  nude  men  in  the 
lateral  and  ventral  positions.  A  value  of  0.0336  ft2 /lb  was  obtained.  This  value  is  higher  than 
the  0.0269  ft2 /lb  obtained  from  the  extrapolation  of  the  animal  data,  but  this  can  be  explained  if 
the  fact  is  taken  into  consideration  that  the  animals  tested  were  much  more  spherically  shaped 
than  man  and  thus  have  a  lower  a  for  the  maximum  presented -area  position. 

Another  way  to  estimate  the  a  of  man  is  by  the  concept  of  the  total  surface  area.  For  a 
sphere  the  total  surface  area  is  47rr2  and  the  average  presented  area  is  7rr2.  For  a  cube  the 
total  surface  area  is  six  times  the  area  of  one  face,  and  the  average  presented  area  of  a  cube 
in  random  orientation  is  %  times  the  area  of  one  face.7  The  ratio  of  the  average  presented 
area  to  the  total  surface  area  is  (7rr2/47rr2)  sphere  =  V4  and  [(3/2)/6]  cube  =  V4* 

In  general,  if  an  object  has  no  cavities  (ellipsoid,  cube,  etc.),  its  average  presented  area 
in  random  orientation  is  one -fourth  the  total  surface  area  (see  Appendix,  Sec.  A. 3,  for  proof). 


Since  a  man  does  have  concavities,  it  is  not  actually  the  total  surface  area  that  should  be  con¬ 
sidered,  but  rather  an  effective  surface  area. 

Bohnenkamp’s  work  with  effective  radiating  surface  area  of  an  erect  man  has  shown  this 
area  to  be  from  83  to  85  per  cent  of  the  total  surface  area.5  The  total  surface  area  can  be  ob¬ 
tained  from  various  sources,8-9  but  for  this  study  the  figures  given  in  the  Schmitt  report6  were 
used  (an  average  of  20.7525  ft2  total  surface  area  and  167.5  lb  weight).  With  these  data  the 
average  value  of  s/m  was  found  to  be  %  x  84%  x  (20.7525  ft2/167.5  lb)  =  0.0260  ft2/lb.  If  C^  is 
assumed  to  be  equal  to  unity,  this  figure  (0.0260  ft2/lb)  can  be  compared  to  the  extrapolated 
value  for  a  of  0.0269  ft2 /lb  (Fig.  4.24). 

The  a  values  derived  from  the  Schmitt  report  for  various  orientations  of  the  men  relative 
to  the  wind  were  used  to  calculate  the  average  a  for  a  rigid,  clothed  man,  giving  a  value  of 
0.0281  ft2 /lb  (see  Appendix,  Sec.  A.4).  The  average  drag  coefficient  for  a  rigid,  clothed  man 
would  then  be  the  average  a  divided  by  the  average  s/m,  or  0.0281/0.0260  =  1.08.  This  value 
compares  very  well  with  the  approximate  values  listed  by  Hoerner,10  which  include  drag  coeffi¬ 
cients  of  1.0  for  a  parachutist  in  free  fall,  1.0  for  a  ski  runner  on  a  slope,  and  1.2  for  a  ski 
jumper  in  air.  It  should  be  remembered  that  the  average  a  of  man  would  be  lowered  if  those 
positions  where  a  man  is  doubled  up,  instead  of  being  rigid,  were  also  considered. 
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TABLE  4.1 — ALPHA  VARIATION  WITH  MASS  DERIVED  FROM  EXPERIMENTAL 

DATA  FOR  SPHERES  AND  STONES* 
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Fig.  4.2  Measured  time  of  fall  minus  no-drag  time  of  fall  for  spheres  vs.  cc  based  on  Cj^  —  0.47.  The  line  is  the  theoretical  curve,  as 
suming  constant  cc  throughout  a  fall  computed  from  Eq.  2.14. 


Plates 
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Fig.  4.3 — Drag  coefficients  vs.  Reynolds  number  for  spheres  and  plates  (maximum  area  presented  to  the  wind).  Included  on  the  figure 
are  the  terminal  points  of  some  of  the  plates  and  spheres  dropped.  The  dashed  lines  represent  average  drag  coefficients  often  used  for  these 
objects  (see  Ref.  3). 
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Fig.  4.4  Effective  <X  vs.  sphere  mass.  The  least -squares  regression  lines  were  computed  from  test  data  assuming  constant  Cp). 
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conditions  as  stated  on  figure. 
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Fig.  4.6  Experimentally  determined  a  vs.  mass  for  limestone  fragments.  Results  of  the  least-squares  analysis  are  shown. 
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from  10P4  and  lSla  locations  at  NTS  are  shown. 


Fig.  4.8  Experimentally  determined  a  vs.  mass  for  IS  la  natural  stones.  Included  are  results  of  the  least-squares  regression  analysis  for 
natural  stones  from  10P4  and  IS  la  locations  at  the  NTS. 


.73833  -  .30 
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Fig.  4.9 — Experimentally  determined  a  vs.  mass  for  natural  stones  from  NTS  at  4.3GTS  location. 
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for  military  debris  showing  results  of  least-squares  analysis.  The  theoretical  line  for 


one  face  down)  computed  from  their  times  of  fall. 


INCHES 


Fig.  4.16  Front  views  of  some  of  the  larger  pieces  of  window  glass  (0.125  in.  thick)  used  in  drop-tests. 
Numbers  indicate  masses  in  grams. 
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Fig.  4.17 — Grouped  data  for  window  glass  dropped  flat  showing  means  and  standard  deviations  and  least -squares  regression  lines  describing 
oc  vs.  mass.  This  figure  indicates  the  need  for  four  straight  line  segments  in  the  analysis. 
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Edge  Window 


Flat  Plate 


Fig.  4.20  —Summary  of  the  glass -fragment  data.  This  shows  regression  lines  describing  ot  vs.  mass  for 
(1)  all  the  glass,  (2)  window  glass  dropped  flat,  (3)  window  glass  dropped  on  edge,  and  (4)  plate  glass 
dropped  flat. 
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RABBIT  1829  g 


Fig.  4.22— “Silhouette  photographs  of  some  of  the  animals  used  in  the  drop  tests.  The  areas  of  the  animals  in  the 
lateral  and  ventral  positions  were  determined  from  full-scale  photographs  of  this  type. 
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Chapter  5 
SUMMARY 


The  purpose  of  this  study  was  to  determine  certain  aerodynamic -drag  parameters  of  small 
irregular  objects,  such  as  window-glass  fragments,  stones,  steel  fragments,  and  small  labora¬ 
tory  animals.  The  drag  parameters  of  interest  were  those  necessary  for  the  computation  of  the 
velocity  vs.  time  history  of  the  object  when  exposed  to  the  winds  associated  with  nuclear - 
produced  blast  waves.1  These  parameters  were  combined  into  one  quantity  called  the  accelera¬ 
tion  coefficient  (a)  and  defined  as  the  product  of  the  area  presented  to  the  wind  and  the  drag 
coefficient  divided  by  the  mass. 

The  experimental  method  used  to  determine  a  was  to  measure  the  time  required  for  a  test 
object  to  free  fall  a  given  distance.  To  relate  a  to  drop  time,  it  was  necessary  to  assume  that, 
first,  the  drag  coefficient  and,  second,  the  area  normal  to  the  object’s  motion  were  both  con¬ 
stant  throughout  the  fall;  i.e.,  a  was  constant.  Experimentally,  the  second  condition  was  not 
difficult  to  satisfy  since  this  only  required  that  the  object  not  rotate.  It  is  well  known  that  the 
drag  coefficient  varies  with  wind  speed;  however,  it  was  shown  from  the  experimental  results 
that  the  error  introduced  in  the  measurement  of  a  by  assuming  that  the  drag  coefficient  was 
constant  was  not  significant,  at  least  in  so  far  as  irregular  objects  were  concerned.  The  follow¬ 
ing  derived  relation  was  used  to  determine  a  from  measured  drop  time: 

T  =  (2H/g/2  (paH)“/2  cosh-1  ePan/2 
where  T  =  drop  time 

g  =  acceleration  of  gravity  (assumed  constant) 

H  =  height  of  drop 

p  -  air  density  (assumed  constant) 

a  =  acceleration  coefficient  (assumed  constant) 

The  equipment  used  in  the  experiment  consisted  of  a  mechanism  that  would  simultaneously 
drop  the  test  object  and  start  an  electronic  timer,  an  impact  plate,  and  a  contact  microphone 
attached  to  the  lower  side  of  the  impact  plate.  When  the  test  object  made  contact  with  the  plate, 
an  electrical  signal  from  the  contact  microphone  stopped  the  timer. 

The  dropping  technique  as  used  in  these  experiments,  was  explored  by  dropping  several 
types  of  regular  objects  (spheres,  cubes,  and  plates);  the  results  obtained  were  compared  with 
those  obtained  using  other  techniques.  The  results  are  presented  in  this  report  for  irregular 
objects.  The  drop  times  for  two  types  of  spheres  were  also  used  to  check  the  experimental 
equipment  for  consistency  of  operation  at  various  times  during  the  drop  tests. 

Several  types  of  stones  were  used  as  test  objects.  Most  of  the  stones  had  been  used  at  the 
NTS  to  study  the  production  of  secondary  missiles  by  blast  winds  produced  by  nuclear  explo¬ 
sions.2,3  From  29  to  45  stones  of  each  type  were  dropped,  and  a  was  determined  for  each.  The 
results  were  analyzed  by  regression  analyses  to  determine  a  for  each  type  of  stone,  and  the 
results  are  listed  below.  The  units  for  a  are  ft2 /lb  and  for  mass  (m),  mg.  Sigma  (a)  is  used  to 
designate  the  standard  error  of  estimate  in  log  units. 
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Gravel  (painted): 


log  a  =  0.57610  -  0.38803  log  m 
30  <  m  <  1500  a  =  0.07463 

Limestone  fragments: 

log  a  =  0.49926  -  0.32159  log  m 
130  <  m  <  15,000  a  =  0.08811 

Natural  stones3  from  NTS  at  location  10P4: 

log  a  =  0.48759  -  0.32835  log  m 
11  <  m  <  4100  a  =  0.07063 

Natural  stones3  from  NTS  at  location  ISla: 

log  a  =  0.45171  -  0.31263  log  m 

17  <  m  <  7100  a  =  0.08684 

Natural  stones3  from  NTS  at  location  4.3  GTS: 

log  a  =  0.73833  -  0.30758  log  m 
11  <  m  <  2300  a  =  0.1294 

Pumice  fragments: 

log  a  =  0.63748  -  0.25466  log  m 

18  <  m  <  2300  a  =  0.08948 

Fragments  of  steel  (mostly  shrapnel)  were  also  used  as  test  objects,  and  the  results  were 
analyzed  in  a  manner  similar  to  that  described  for  stones.  The  resulting  regression  equation 
was 

log  cl  =-0.74223  -  0.11315  log  (m/1000) 

4500  <  m  <  150,000  a  =  0.09974 

Fragments  of  ordinary  window  glass  0.125  in.  thick  were  dropped,  both  with  a  plane  sur¬ 
face  down  and  edge -on,  to  determine  the  maximum  and  minimum  a  as  a  function  of  mass.  In 
addition,  plate -glass  fragments  0.225  in.  thick  were  dropped  with  the  plane  surface  down.  For 
small  fragments  the  a  values  determined  were  independent  of  both  orientation  and  type  of  glass 
(ordinary  window  glass  or  plate  glass).  It  was  found  that  the  relation  between  a  and  mass 
could  be  represented  as  a  straight  line  on  logarithmic  paper,  but  only  for  definite  mass  in¬ 
tervals.  The  regression  equations  derived  from  the  experimental  data  are  given  below  (n 
represents  the  number  of  fragments  dropped;  the  remaining  nomenclature  is  the  same  as  that 
used  previously). 

Window  glass  and  plate  glass  dropped  flat  or  edgewise  for  10  <  m  <  220: 

log  a  =  0.55659  -  (1/3)  log  m 
n  =  155  ct  =  0.09014 
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Window  glass  dropped  flat  for  220  <  m  <  1270: 

log  a  — -0.13874  -  0.036374  log  m 
n  =  81  a  =0.06348 

Window  glass  dropped  flat  for  1270  <  m  <  1574: 

log  a  =  -3.02323  +  0.89295  log  m 
n  =  10  a  =  0.04525 

Window  glass  dropped  flat  for  1574  <  m  <  105: 

log  a  =  -0.28024  +  0.34985  log  m 
n  =  46  a  =0.01473 

Window  glass  dropped  edgewise  for  220  <  m  <  5000: 

log  a  =  0.12689  -  0.14982  log  m 
n  =  88  a  =  0.1013 

Plate  glass  dropped  flat  for  860  <  m  <  3  x  105: 

log  a  =  0.49172  +  0.023921  log  m 
n  =  31  a  =  0.03416 

A  useful  formula  was  derived  which  related  drag  coefficient  (Cd)  for  plates  of  random 
shapes  and  finite  thickness  to  the  ratio  of  the  square  root  of  plate  surface  area  (s)  to  thickness 
(t).  The  drop-test  data  used  in  this  analysis  included  that  for  window-  and  plate-glass  frag¬ 
ments  dropped  flat  as  well  as  that  for  plates  of  regular  shapes  (square,  circular,  triangular) 
made  of  brass  and  plastic.  The  regression  equation  noted  below  was  computed  from  data  for 
117  test  objects  whose  shape  factor  (s1/2/t)  varied  from  2  to  65. 

log  CD  =  0.019445  +  0.048498  log  (Sl/2/t) 
a  =  0.02428 

(Note:  s1/2  and  t  are  expressed  in  the  same  units.) 

Shadowgraphs  were  prepared  for  mice,  rats,  guinea  pigs,  and  rabbits  in  ventral  and  lateral 
positions  to  determine  their  presenting  areas.  In  general,  the  areas  determined  (Fig.  4.23) 
were  proportional  to  the  mass  of  animal  raised  to  the  2/3  power.  The  animals  were  then 
dropped  ventrally  and  laterally  to  determine  their  a.  It  was  found  that  there  was  no  significant 
difference  in  a  determined  for  these  positions.  The  experimentally  determined  a  values  for 
each  animal  species  were: 

Mice  (20.44  g)  0.3769  ±  0.0425  ft2/lb 

Rats  (184.48  g)  0.1853  ±  0.0142  ft2 /lb 

Guinea  pigs  (530.1  g)  0.1547  ±  0.0196  ft2/lb 
Rabbits  (2087  g)  0.0794  ±  0.0155  ft2 /lb 

These  data  were  used  to  determine  a  as  a  function  of  animal  mass. 

log  a  =  0.01153  -  0.32400  log  (m/1000) 
n  =  84  or  =  0.06496 
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An  estimate  of  the  a  for  a  168-lb  man  was  determined  to  be  0.0269  ft2/lb  by  the  above 
equation.  Another  estimate  was  made  for  the  a  of  a  16 8 -lb  man  using  the  results  of  a  wind- 
tunnel  study  reported  by  Schmitt.4  For  positions  corresponding  to  those  of  the  dropped  animals 
(ventral  and  lateral),  an  average  a  of  0.0336  ft2 /lb  was  obtained. 

A  method  for  estimating  the  average  presenting  area  of  any  object  in  random  orientations 
was  developed.  It  was  shown  that  the  average  presenting  area  is  one-fourth  the  total  surface 
area,  provided  no  cavities  exist.  To  approximate  the  no-cavity  condition,  the  radiating  surface 
area  for  man  was  used  to  estimate  the  average  a  for  the  human  subjects  used  in  the  wind- 
tunnel  study  mentioned  above.  Assuming  a  drag  coefficient  of  1.0,  a  value  of  0.0260  ft2 /lb  for  a 
was  obtained.  This  value  compared  well  with  a  value  of  0.0281  ft2/lb  obtained  by  averaging  the 
wind-tunnel  data  for  men  in  various  positions  relative  to  the  wind. 
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Appendix 

DERIVATION  OF  EQUATIONS 


A.l  EQUATIONS  OF  MOTION 

A.  1.1  Integration  of  the  Basic  Equation 
Equation  2.9  can  be  integrated 
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Now,  solving  Eq.  A. 2  for  v2  gives 
J/2 
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(A.  3) 


where  h  =  h2  -  represents  the  height  fallen  from  the  point  where  the  velocity  equalled  v*. 
Substituting  v2  =  dh/dt  in  Eq.  A.  3  and  integrating  from  zero  height  of  drop  to  h 
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gives 

Substituting  h  =  (h2  -  ht)  from  Eq.  A.2  into  Eq.  A.5  gives 
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A.  1.2  Dimensionless  Analysis 


Equations  2.4  and  2.8  of  the  text  can  be  written  in  the  form 


dv  1  2  s  p/vD  \ 

v  — -g--pv2  — f(  —  \ 


dh 


m  V  v  ) 


(A.7) 


where  f(vD/y)  =  f(Rd)  =  Cd,  the  drag  coefficient. 

Now,  g,  D,  and  v  were  considered  constant  so  that  they  could  be  placed  inside  the  differ¬ 
entials.  Algebraic  manipulations  were  then  made  on  Eq.  A.7  to  make  dimensionless  quantities 
as  follows: 


vD 


d(-) 
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=  1 


Defining  Y  =  hg(D/V)2  (distance  numeric) 

Z  =  vD/v  (velocity  numeric) 

K  ^  (p/g)  (s/m)  (^/D)2  (acceleration- coefficient  numeric) 

as  numeric  parameters,  the  above  dimensionless  equation  becomes 


(A.  7a) 


zdf  =  1-^Kz2  f(z) 


(A. 8) 


By  basic  definition,  V  =  AH/AT  or  AT  =  AH/V,  Eq.  A. 8  does  not  contain  the  time  numeric. 
However,  the  incremental  step  of  time  numeric,  defined  as  AX,  would  be  the  incremental  step 
of  distance  numeric  divided  by  the  velocity  number.  Restated  in  terms  of  an  equation,  this  is 
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(A. 9) 


Any  previously  derived  equation  can  now  be  converted  to  the  nondimensional  form  by  re¬ 
placing  the  primary  parameters  v,  h,  t,  and  s/m,  by  the  numerics  Z,  Y,  X,  and  K.  The  p,  v , 
and  g  quantities  are  dropped  in  the  nondimensional  form  as  they  go  into  the  numerics.  CD  re¬ 
mains  the  same  in  the  nondimensional  equations. 

Equations  A. 2  and  A. 6  become 
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(T  -  T^d)  or  time-of-fall  minus  no-drag  time-of-fall  parameter*  (Sec.  2.5),  would  now  be 


T  -  Tnd  =  T  - 


X  —  (2Y) 


Another  useful  parameter  (Sec.  2.5)  would  be 


(A. 12) 


paH 

JT  AND 


T  —  T  d  —  ^  ~  T  vn  —  T 


2H  \1/2 


/2H 

\T 


1  + 


P-CD  H 
m 

12 


-  X  -  (2Y) 


(>+%ky) 


(A. 13) 


If  the  approximate  average  value  of  the  drag  coefficient  of  an  object  is  substituted  for 
in  Eq.  A.  13,  the  parameter  then  represents  the  scaled  time  difference  between  the  actual  time 
above  no -drag  time  and  an  approximation  to  the  time  above  no-drag  time  considering  the  drag 
coefficient  to  be  constant  at  Cd-  It  should  be  noted  one  can  go  from  the  variable  of  Eq.  A. 13  to 
that  of  Eq.  A.  12  by  adding  2Y(Cd/12)  KY,  which  is  the  approximation  to  the  time  of  fall  minus 

no-drag  time  of  fall  in  nondimensional  form. 

The  no-drag  equations,  Eqs.  2.18  and  2.19,  can  likewise  be  put  in  numeric  form  to  give 
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Equations  A.  10  and  A.  11  are  the  equations  actually  used  to  predict  the  accurate  drop  times 
of  two  spheres  (Sec.  4.2).  The  CD  curve  was  broken  up  into  about  one  hundred  steps  over  the 
Reynolds -number  range  encountered  by  these  spheres,  and  the  CD  values  were  used  to  integrate 
step  by  step  to  the  total  height  of  drop. 


A. 2  TERMINAL  VELOCITY 


Although  not  apropos  of  the  present  experiment,  it  should  be  noted  that  a  useful  relation 
can  be  derived  from  Eq.  2.13.  As  a  body  falls  farther  and  farther,  its  velocity  tends  to  reach  a 
finite  value  called  terminal  velocity.  This  terminal  velocity  can  be  found  from  Eq.  2.13  by 
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Figure  A.l  is  a  plot  of  terminal  velocity  vs.  a  for  two  different  conditions  of  p  and  g. 

The  percentage  of  terminal  velocity  reached  by  an  object  with  constant  a  after  it  has  fallen 
a  distance  H  can  be  computed  from  Eqs.  3.13  and  A.  16  by 
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Figure  A. 2  is  a  plot  of  the  above  relation  showing  the  percentage  of  terminal  velocity  as  a 
function  of  a  and  H.  Equation  A.  16  is  independent  of  any  initial  conditions  and  represents  the 


♦The  subscript  ND  is  used  to  indicate  the  no-drag  value  of  a  parameter;  i.e.,  TND  is  the 
time  an  object  would  take  to  fall,  assuming  the  object  did  not  experience  drag  force. 
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velocity  to  which  the  instantaneous  velocity  is  tending  at  any  given  time;  whereas,  Eq.  A. 17 
requires  that  the  velocity  be  zero  when  the  height  is  zero  and  also,  that  p,  a,  and  g  remain 
constant  throughout  the  fall. 

A. 3  AVERAGE  PRESENTED  AREA 

Consider  a  smooth  surface,  S,  which  has  no  cavities  (part  a,  Fig.  A. 3);  i.e.,  a  plane  that  is 
tangent  to  S  at  any  point  does  not  cut  through  the  surface.  Then  at  some  point  on  S  construct  a 
tangent  plane  that  approximates  S  over  a  small  area  dA  in  the  vicinity  of  the  point  of  tangency. 
Let  N  be  the  normal  to  the  plane  at  the  point  of  tangency,  E  be  oriented  in  any  reference  direc¬ 
tion  fixed  in  space,  and  9  be  the  angle  between  N  and  E.  As  N  takes  random  orientations,  not 
all  values  of  9  are  equally  likely.  To  find  the  probability  distribution  of  9 ,  one  compares  the 
area  on  the  surface  of  a  sphere  between  9  and  9  +  d0  to  the  total  area  of  the  sphere.  That  is  to 
say,  consider  a  sphere  whose  center  is  always  located  at  the  point  of  tangency  of  the  plane  on 
the  smooth  surfaces.  Now  N  can  intersect  the  sphere  in  any  point  on  the  sphere,  and  all  such 
points  are  equally  likely.  Then  the  probability  that  N  will  intersect  the  sphere  in  a  specified 
area  of  the  sphere  is  just  the  ratio  of  this  area  to  the  total  area  of  the  sphere.  In  the  case 
under  consideration,  the  specified  area  is  the  total  area  on  the  sphere  between  the  angular 
limits  of  9  and  9  +'d0  with  respect  to  the  fixed  reference  direction  E.  This  area  (see  part  b, 
Fig.  A. 3)  is  27r(r  sin  9)  r  69;  so  the  probability  that  the  angle  between  N  and  E  is  between  9  and 
9  +  d0  is 


m  de  =  ^(TSlnJ)rde 


(A. 18) 


The  presented  area  in  the  direction  of  E  produced  by  dA  (part  a,  Fig.  A. 3)  is 


Presented  area  =  dA  cos  9 
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For  7t/2  <  9  <  tt,  dA  could  not  be  seen  from  the  direction  of  E  because  it  is  behind  S.  The  aver¬ 
age  presented  area  would  then  be  the  integral  of  the  presented  area  as  a  function  of  9  weighted 
with  the  probability  distribution  of  9,  or 
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For  the  whole  surface  S,  the  average  presented  area  divided  by  the  total  surface  area  would  be 


Average  presented  area 
Total  surface  area 


fsdA_i 

I 


(A. 21) 


A. 4  AVERAGE  a  FOR  A  RIGID  MAN 


Two  angles  were  used  to  fix  the  man  with  respect  to  the  wind,  9  and  <p  (part  c,  Fig.  A. 3). 
Alpha  can  be  computed  directly  from  Schmitt’s  report1  for  9  =  90°;  0°  <  (p  <  360°  (upright 
position)  or  for  0°  <  9  <  180°;  <p  =  90°  (side  to  the  wind).  For  intermediate  values  the  following 
formula  was  used  to  approximate  ol 


a(0,<t>)  «  «(0°)  +  [«<90»  -  g(0°)J  [«<*»  90D>-«<0°>] 

a(90°,90°)  -  a(0°) 


(A. 22) 
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where  a(0°)  represents  the  a  of  a  man  head-on  to  the  wind  (9  —  0)  and  is  the  same  irrespective 
of  <p,  similarly  a(180°)  represents  the  a  of  a  man  feet-on  to  the  wind  (9  =  180°)  and  also  is 
constant  for  all  values  of  <p.  From  Schmitt’s  report1  6(0°)  =  6(180°)  for  a  clothed  standing  man 
and  has  a  value  of  0.00627  ft2 /lb.  Thus  Eq.  A.22  holds  exactly  for  all  the  conditions  that  can  be 
computed  directly  from  Schmitt’s  report;  i.e.,  substituting  in  Eq.  A.22  the  various  values  of  9 
and  <p  for  which  a  is  known,  the  following  relations  for  a  are  obtained: 


a(0  °,cp)  =  a(0°) 
a(180  °,cp)  =  a(0°) 
a(90°,<p )  =  a(90  °,cp) 
a(6,90°)  =  a(@,90°) 


(A. 23) 


Equation  A.22  was  arrived  at  by  assuming  [a(6,c/>i)  -  a(O)]/[a(0,02)  -  a(0)]  is  constant  for  all 
values  of  9  when  cpt  and  <fo  are  fixed  or,  alternatively,  by  assuming  [a(9u<p)  -  a(0°)]/[«(62,<#>)  - 
a(0°)]  is  constant  for  all  values  of  (p  when  0,  and  02  are  fixed.  In  connection  with  the  first 
assumption,  one  can  write 


a(e,<p)-a( 0°)  _  o>(9O°,0)  -  ot(0°)  (A  24) 

a(0,9O°)  -  a(0°)  a(90°,90°)  -  a( 0°) 

for  all  values  of  6  and  cp.  Solving  this  for  a:(0,(f>)  gives  Eq.  A.22.  The  second  assumption  will 
yield  Eq.  A.22  in  a  similar  manner.  Since  6  and  (p  are  independent  variables,  the  average  a 
would  be  given  by 


a{Q,(p)  =  a^O0)  + 


[a(90 °,(p)  -  a(0°)l[a(S,90o)  -  ct( 0°)] 
a(90°,90°)  -  a(0^) 


(A. 25) 


where  the  bars  indicate  average  values. 

The  value  a(0,9O°)  was  computed  by  weighing  a(0,9O°)  with  the  probability  distribution 
(sin  6/2)  d0  (Eq.  A. 3)  and  integrating  graphically  from  6  =  0  to  6  =  n  radians.  The  value  for 
o/(90°, <p j  is  just  the  average  value  of  a(90°,^j  with  no  weighing  (all  values  of  <p  are  equally 
likely),  and  a(90°,90°)  and  a(0°)  are  particular  a  values  computed  directly  from  Schmitt's  re¬ 
port.  The  value  for  a(6,(p)  computed  from  Eq.  A. 25  represents  the  average  a  for  a  rigid, 
straight  man,  assuming  every  possible  orientation  with  respect  to  the  wind  is  equally  probable. 
(Average  a  computed  from  Schmitt's  report  is  0.0281  ft2 /lb.) 
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Terminal  velocity  vs.  a  for  standard  conditions  at  sea  level  (curve  I)  and  typical  conditions  at  5300-ft  altitude  (curve  II). 
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Fig.  A.2— Percentage  of  terminal  velocity  vs.  height  of  fall  and  a.  When  p  =  0.0653  lb/ft3  (typical  air  density  for  5300  ft  above  sea  level), 
the  height  can  be  read  directly,  if  p  does  not  have  this  valuej  the  height  variable  must  be  multiplied  by  the  appropriate  factor  to  convert  to 
actual  height  of  fall. 
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Through  its  Division  of  Biology  and  Medicine  and  Civil  Effects  Test  Opera¬ 
tions  Office,  the  Atomic  Energy  Commission  conducts  certain  technical  tests, 
exercises,  surveys,  and  research  directed  primarily  toward  practical  applica¬ 
tions  of  nuclear  effects  information  and  toward  encouraging  better  technical, 
professional,  and  public  understanding  and  utilization  of  the  vast  body  of  facts 
useful  in  the  design  of  countermeasures  against  weapons  effects.  The  activities 
carried  out  in  these  studies  do  not  require  nuclear  detonations. 

A  complete  listing  of  all  the  studies  now  underway  is  impossible  in  the 
space  available  here.  However,  the  following  is  a  list  of  all  reports  available 
from  studies  that  have  been  completed.  All  reports  listed  are  available  from 
the  Office  of  Technical  Services,  Department  of  Commerce,  Washington  25, 

D.  C.,  at  the  prices  indicated. 
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